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ABSTRACT 

We present the first results of the SOAR Gravitational Arc Survey (SOGRAS). The 
survey imaged 51 clusters in two narrow redshift intervals centered at z — 0.27 and z = 0.55, 
targeting the richest clusters in each interval. Images were obtained in the g', r', and i' bands 
using the SOAR Optical Imager (SOI), with a median seeing of 0.82", 0.74", and 0.69", 
respectively, in these filters. Most of the survey clusters are located within the Sloan Digital 
Sky Survey (SDSS) Stripe 82 region and all of them are in the SDSS footprint. Photometric 
calibration was therefore performed using SDSS stars located in our SOI fields. We reached 
for galaxies in all fields the detection limits of g ~ 23, r ~ 22.5, and i ~ 22 for S/N = 3. 
As a by-product of the image processing, we generated a source catalog with 21280 entries, 
the vast majority of which are galaxies, where we list their positions, magnitudes, and shape 
parameters. We compared our galaxy shape measurements to those of local galaxies and 
concluded that they were not strongly affected by seeing. From the catalog data, we are able 
to identify a red sequence of galaxies in most clusters in the lower z range. We found 16 
gravitational arc candidates around 8 clusters in our sample. They tend to be bluer than the 
central galaxies in the lensing cluster. A preliminary analysis indicates that ~ 10% of the 
clusters have arcs around them, with a possible indication of a larger efficiency associated to 
the high-z systems when compared to the low-z ones. Deeper follow-up images with Gemini 
strengthen the case for the strong lensing nature of the candidates found in this survey. 
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1 INTRODUCTION 

Tracing the evolution of galaxy cluster properties, in particular their 
mass distribution, has important implications for their use as cos- 
mological probes, for understanding the nature of dark matter and 
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dark energy, and to constrain galaxy evolution. A unique way to 
assess the mass distribution in clusters is through the arcs produced 
by strong gravitational lensing (Blandf ord & Narayan|1992{|Hat-| 
|tori et al.|[T999l |Guzik & Seljak||2002[ |Mandelbaum et al.||2006| 
|Treu| |2010; Kneib & Natarajan 2012). In particular, the statistics 
of gravitational arcs may provide constraints on cosmological pa- 
rameters and on scenarios of structure formation (Bartelmann et al. 
|1998||2003l|Golse, Kneib & Soucail|2002l|Meneghetti et al.|2004l 



IKochanek etalT2 006 Hilb ert, Metcalf & W hite 2007| |Vuissozet| 
|al.|2007) . 

This motivated arc searches to be conducted, both in images 
from wide field surveys (Gladders et al.|2003||Estrada et al.|2007] 
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ICabanac et a l. 2007 ; Beloku roTl* al.|2009t|Kubo et al.|2010||Kneibl 
[eTaT1[20T0l |Gilbank et aT1[20TT] |More et al.|20l2"l >, as well as in 
fields targeting known clusters JLuppino et al.||1999| |Zaritsky &| 
|Gonzalez|[2003l |Hennawi et al.[|2008| |Kausch et al.||2010[ >. Up- 
coming wi de field imaging surveys, such as the Dark Energy Sur 



H(DES; 



Annis et al. 



2005 



The Dark Energy Survey Collabora- 



|tion|2005) , which starts operations this year, are expected to detect 
about an order of magnitude more arcs than the current largest sur- 
veys. 

For over a decade there has been a debate about the compat- 
ibility of the observed arc abundance with theoretical modeling. 
|Bartelmann et al.| ( |1998fr suggested an apparent overabundance by 
approximately an order of magnitude of giant arcs on the sky as 
compared to ACDM predictions. Subsequent comparison works 
with limited statistics confirmed that the number of giant arcs on 
the sky is underpredicted by the ACDM cosmological model ( |Lup-| 
pino et al.|f999l|Gladders et al.|2003"||Zaritsky & Gonzalez|2003| 
Li et al.|200oy 

However, more recent studies including several factors that 
were not present in the first predictions and looking at the appro- 
priate observable (e.g., the number of arcs per cluster) have reduced 
the overall discrepancy between observed and predicted arc abun- 
dances ( [Dalai, Holder & Hennawi|2004| [Horesh et al.|2005l[Herri 
|nawi et al.|2 007 ). In particular, Horesh et al. (20lTjl have carried 
out a study using simulations at the image level, where gravita- 
tional arcs are generated from ray-tracing of realistic sources from 
the Hubble Ultra Deep Field through clusters from the Millennium 
A-body simulation. The fraction of arcs per cluster in the simu- 
lated samples is compared to an arc sample in clusters with simi- 
lar properties, using the same methods to identify the arcs in both 
samples. They find an overall consistency of the observed and sim- 
ulated samples, at least in the redshift interval 0.3 < z < 0.6. 

If on the one hand the arcs statistics problem may have been 
solved or at least mitigated, issues remain regarding the variation 
of arc abundance with respect to the cluster redshift. For example, 
Gladd ers et aT] p003] > found an over-abundance of arcs in high- 
redshift clusters as compared to lower redshift ones. Gon zalez et| 
[aLl pOT2l found arcs in a cluster at z = 1.75, which should not be 
present at their image depths according to their modeling. Horesh 
|et al.|p011} found an under-prediction of clusters at z ~ 0.2 as 
compared to observations analyzed in Horesh et al. (20051. These 
discrepancies could be due to the evolution of cluster structure with 
redshift (including the role of baryons) and/or to selection effects 
of the samples. Cam inha et al.| (|2012) model the variation of arc 
abundance with lens redshift stressing the effect of magnification 
on the expected distribution and finding an increase of arc incidence 
with z. 

The main motivation for the SOAR Gravitational Arc Survey 
(SOGRAS) is to constrain the variation of strong lensing efficiency 
as a function of cluster redshift, comparing the results with theo- 
retical expectations. For this sake, we have designed a survey tar- 
geting clusters distributed in two redshift bins centered at z ~ 0.3 
and z ~ 0.5. A total of 51 clusters were imaged in the g' , r', and 
i' bands with the 4. 1 Southern Astrophysical Research Telescope 
(SOAR) from mid 2008 to the end of 2010. 

The arcs and other strong lensing features found can be used 
to constrain the individual masses of the clusters (e.g., |Cypriano| 
|et al.|2005) . Another valuable information that can be drawn from 
the data is an estimate of ensemble cluster masses in each z bin 



with weak lensing, by stacking the profile of the scaled tangential 
distortion of background sources of all clusters in that bin. This 
technique has been applied successfully ( Sheldon et al. 200 1 j |2004[ 
|Johnston et al. 2007) and leads to an averaged overall mass for the 
clusters. 

Another motivation was to use this dataset as a test bed for 
tools being developed for DES, in particular for gravitational arc 
studies, including testing automated arc-finders (using either mor- 
phology or colour) and methods to measure arcs properties (e.g. 
|Furlanetto et al.|20 12). Indeed SOGRAS has comparable depth and 
seeing conditions as expected from DES and covers 3 of the 5 DES 
bands. 

Finally, SOGRAS can be seen as a pathfinder for a high resolu- 
tion arc survey with S OAR using the recently commissioned SOAR 
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Adaptive Modultp]( Tokovinin & Cantarutti|2008 
|20T0) , 

In this paper we report on the overall properties of the survey, 
from the target selection and observations to data reduction and 
photometric calibration. We present the photometric catalogs, and 
discuss new gravitational arc candidates found by visual inspection. 
Detailed results on arc analyses and comparison with theoretical 
modeling will be presented in a accompanying paper. 

As a by-product of the survey a large catalog of galaxies in 
the cluster fields was generated, with astrometric, photometric, and 
morphological information. This catalog was used to separate clus- 
ter members (through the red sequence in the colour-magnitude di- 
agrams) from field galaxies, and will be useful for future analysis 
of galaxy evolution. 

The outline of this paper is as follows: ^describes the sur- 
vey, including information on sample selection, observational and 
image details. In [|3]we describe the data reduction, including as- 
trometric and photometric calibrations. We also carefully assess the 
quality of our photometry. The resulting galaxy catalog is presented 
in ij4] In the same section we present the first sample of arc candi- 
dates. Finally, in 5j5]we present our summary and closing remarks. 



2 THE SURVEY 

We have designed the survey to image a sample of galaxy clus- 
ters, equally split into two disconnected redshift bins, one at 
0.20 < Zphot < 0.35 (the "low-z" bin) and the other at 
0.50 < z phot < 0.60 (the "high-z" bin), to have a "leverage 
arm" to constrain the evolution of arc incidence between these two 
intervals. 

The low-z bin was chosen such that there are enough reason- 
ably rich clusters in this bin on the survey footprint (see |2.1[ l and 
to avoid having a too small arc probability. The high-z bin was 
determined by the availability of optical cluster catalogs in the sur- 
vey footprint and by the requirement of having enough background 
sources to allow for a weak lensing analysis by stacking the clusters 
in this bin (to have an overall estimate of the cluster masses). 

All clusters were observed with the SOAR telescope, located 
on Cerro Pachon in the Chilean Andes, with the SOAR Optical 
Imager (SOI). The choice of telescope and instrument is motivated 
by the typical site seeing (~ 0.8") and detector pixel size, which 
yield the required image quality for gravitational arc detection^] 



www .ctio.noao. edu/new/ Teles copes /SOAR/ Instrument s/ SAM/ 
3 As is well known, the detectability of gravitational arcs is very sensitive 
to the PSF FWHM, because the seeing tends to decrease their length-to- 
width ratios and dilutes their surface brightness (e.g. |Cypriano et al.|2001| . 
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The observations of all our targets were carried out in queue- 
scheduled mode assuring that our quality requirements were met. 
Therefore, this survey provides a fairly homogeneous sample, in the 
sense that all images were obtained in similar conditions, with the 
same instrument and filters and same exposure time, and is there- 
fore well suited for a comparison of arc incidence. 

The exposure time was determined by a balance between the 
number of clusters to be observed and the depth achieved for each 
cluster field for a given total observing time. Using the model for 
the number of arcs expected per cluster (as a function of limiting 
magnitude, cluster redshift, etc.) given in Caminha et al. (2012} and 
the exposure time calculate^] we found that the maximum total 
number of arcs is reached for integration times of about 10 min (in 
one single filter). 

We also required to image in 3 bands such that colour infor- 
mation could be gathered, since this is an important information 
for discriminating gravitational arcs from cluster tidal features and 
for identifying multiple images. Furthermore, the colour informa- 
tion helps mitigating the contamination by foreground objects for 
weak-lensing mass reconstructions. 



2.1 SDSS and Stripe 82 

The baseline footprint for the targeting of our survey was the SDSS 
"Stripe 82", a 275 deg 2 equatorial stripe (over —50 < RA < 59; 
— 1.25 < DEC < 1.25), which was scanned multiple times in the 
fall seasons of 2000-2007 as part of a supernovae search, leading to 
a much deeper survey. The final Stripe 82 coadded data (hereafter 
coadd) reaches r ~ 23.5 for galaxies, i.e. 2 mag fainter than the 
main SDSS survey ( [Annis et al.|201 1( >. 

The availability of these data allowed for the construction of 
deeper cluster catalogs, well suited for our target selection. In par- 
ticular, red-sequence based cluster catalogs started to be produced 
as the first coadds of several stripe 82 visits were created. While sin- 
gle pass SDSS cluster catalogs reached up to z ~ 0.3 ( Koes ter et al.| 
|2007| >, catalogs obtained from the coadd reach z ~ 0.6, matching 
our requirement for the high-z bin, having at the same time some 
leverage arm with respect to the low-z bin and still allowing for a 
stacked weak lensing measurement from our data. 

At present, besides the deeper SDSS imaging, parts of Stripe 
82 have been covered by a wealth of multi-wavelength data, such 
as the UKIDSS Large Area Survey in the YJHK bands (|Lawrence| 
[eTaT1[2007), deep GALEX UV imaging JMartin et al.||2005>, the 
SHELA ( |Papovich et al.|2011| > and SpIES (Richards et al.|20l2) 
surveys with Spitzer/IRAC, and the Herschel/HerMES Large Mode 
Survey {Oliver et al.||2012| >. At longer wavelengths, the whole 
Stripe lies within the footprint of Atac ama Cosmology Telescope 
equatorial survey I Sehgal et al.|2012| and 80 deg 2 of the Stripe 
have deep VLA data ( |Hodge et al" 2011} . Stripe 82 has also a 
very high density of spectroscopic redshifts, with redshift measure- 
ments from SDSS ( |Abazajian et al.|2009|l, 2dF <|Colless et al.|2001| 



Croom etal.|2001), 2SLAQ (|Croom, et al. 2009), 6dF (Jones et al.| 



2009l >, DEEP2 (Newman et al.|2012) , VVDS ( |Garilli et al.|2008> , 



PRIMUS jCoil et al. 201 1 ), SDSS-III/BOSS ( |Ahn et al.|20T2) , and 
WiggleZ 1 Drink water et a [ 20 1 ) . This region is thus emerging as 
a "deep extragalactic survey field", a precursor to DES and LSST, 
with an impressive array of multi-wavelength observations already 
in hand or in progress. SOGRAS may be used as a test case for 
combining the good imaging data from SOAR with this large set of 



complementary data, which strengthens the case to carry out most 
of our selection in this field. 

In particular, during the second semester of 2010 until early 
2011, 170 deg 2 of Stripe 82 were imaged in the i band for the 
CFHT/Megacam Stripe 82 Survey (CS82, |Knei"b et al.|2010||Erben| 
|et al.|2012l >, providing data down to i = 23.5 obtained in excellent 
seeing conditions (median seeing of 0.6"), enabling precision weak 
lensing measurements. This dataset is particularly synergistic with 
SOGRAS. CS82 data will allow us to study the clusters imaged for 
SOGRAS at much larger radii and provide weak lensing measure- 
ments for them. On the other hand the SOGRAS data is useful for 
quality assessment on this new dataset around clusters. For exam- 
ple, the star-galaxy separation and the determination of background 
and foreground sources can be tested in these fields thanks to the 
color information at higher depths and better seeing than the SDSS 
photometric data. 



2.2 Selection of galaxy cluster sample 

The survey was carried out in two seasons, the first during semester 
2008B and the second in 2010B JMakler et al.|2008|[20T0| >. Since 
the cluster catalogs and status of the Stripe 82 coadd evolved during 
the two seasons, different catalogs were used for the selection. The 
procedure was nevertheless the same for both seasons: selecting the 
richest cluster catalogs in the same two redshift bins and requiring 
the same imaging conditions and instrument configurations. There- 
fore the two sets of observations are considered as a single dataset. 

The cluster selection for the 2008B season was carried out us- 
ing a combination of three unpublished cluster catalogs on Stripe 
82 (J. Hao, T. McKay, et al.). The cluster finding methods were 
based on the red-sequence, accounting for its variation with red- 
shift, and are precursors of the Gaussian Mixture Brightest Cluster 
Galaxy (GMBCG) cluster finder ( |Hao et al.|2010| >. However they 
differ in their likelihoods and the radial profiles used. All were run 
in the coadded data available in late 2006, providing an estimate 
of the cluster photometric redshift {z v hot> and richness (roughly 
the number of red-sequence galaxies in the cluster with luminosity 
above L* /2). We selected the richest clusters from these catalogs 
in the two redshift bins and ranked them by richness. 

We complemented this main sample, with an extra sample 
consisting of clusters detected on SDSS Data Release (DR) 6 data 
l |Adelman-McCarthy et al.||2008|l — n ot necessarily on stripe 82 
— from the MaxBCG {Koester et al.||2007} catalog that matched 
ROSAT x-ray sources and had good observability from SOAR on 
that semester. No redshift restriction was applied to this sample, 
which was chosen to improve our chances of finding arc systems 
and for scheduling flexibility (i.e. to allow observations to be made 
when observing conditions were not suitable for Stripe 82). Natu- 
rally, clusters observed from this sample are not suitable for our arc 
statistics purposes. 

A visual inspection of SDSS single pass images using the Cat- 
alog Archive Serve^] was made in order to avoid clusters close 
to bright stars, which could jeopardize the observations. We dis- 
carded all clusters that show diffraction spikes and star halos within 
a ~ 6.5' x 6.5' field around the cluster center. We also discarded 
fields with saturated stars within 3' from the cluster center, im- 
posing stronger limits on the magnitude closer to the center (e.g. 
mag < 14 for 6 < 1'). This eliminated ~ 25% of the selected 
fields. Clusters with significant substructure, or which appeared to 



4 http : / / www .noao.edu/ gateway /ccdtime/ 



5 http://cas.5d5s. org /dr6/en/tools/ chart /list .asp 
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be line-of-sight superpositions or spurious detections from this vi- 
sual inspection were also avoided, eliminating ~ 10% of the se- 
lected clusters. While performing the visual inspection, we ignored 
any potential arc feature to avoid biasing the sample. The final re- 
sult of this process was a set of two lists (one for each redshift bin) 
ordered by richness, containing a total of 60 selected clusters. The 
observers were told to select randomly among these lists, choosing 
the highest ranked object for which the observing conditions were 
favorable. 

A total of 18 clusters were observed in that season, 13 in the 
high-z bin, 4 in the low-z bin, and 1 from the extra sample. 

For the 201 OB season the clusters were selected from a GM- 
BCG cluster catalog constructed using the complete coadd data. To 
exploit the synergy with CS82, only clusters in the footprint of this 
survey were selected. Again we selected the clusters on the two 
redshift bins and kept the richest ones. As in the 2008 sample, we 
also included an extra sample with the same objects selected for 
that season. The visual selection procedure was the same as for the 
2008 season. 

Initially, 26 clusters were observed, 12 in the high-z bin, 1 1 in 
the low-z bin, and 3 from the extra sample. At that point, the SO- 
GRAS program had still telescope time allocated, but the observ- 
ability of Stripe 82 was unfavorable. Since there were more clusters 
observed from the high-z bin, than in the low-z one, clusters could 
be selected only in the later, thus requiring only shallower imaging. 
Therefore an auxiliary sample sample was chosen, following the 
same selection criteria as the low-z, but choosing clusters at higher 
RA in an equatorial region covered by SDSS single pass imaging. 
These clusters were taken from a GMBCG catalog based on SDSS 
DR7 data ( |Hao et aL]|2010| >. 7 clusters were observed from this 
sample, completing the survey. 

A summary of the properties of SOGRAS observed clusters is 
given in Table[T] The on-sky distribution of these clusters is shown 
in Figure[T] From the 51 observed clusters, 40 are in the main sam- 
ple, 7 are in the auxiliary 2010 sample and 4 are in the extra sample. 
For arc statistics analysis, the auxiliary sample can be added to the 
main sample, since their clusters were selected following the same 
criteria for redshift and richness. Although the clusters in the aux- 
iliary sample were taken from a shallower photometric catalog, the 
same cluster finder algorithm (GMBCG) was used to generate the 
cluster catalog in the main SDSS area and in the coadd and it is 
known to be complete for the low-z clusters. The resulting split 
in z was 25 clusters in the high-z bin and 22 in the low-z bin. 
The 4 clusters in the extra sample cannot be included in the statis- 
tical analysis because they were selected following other criteria. 
The distribution of photometric redshifts for the SOGRAS clusters, 
taken from the catalogs used for the target selection, is shown in 
Figure [2] 

2.3 Observations 

Observations of all our targets were carried out with SOI, which 
consists of a mini-mosaic of two E2V CCDs, each one with 
4096 x 2048 pixels, covering a field of view of 5.25' x 5.25'. A 
2x2 binning was used, yielding a detector scale of 0.154"/pixel. 
The exposures were taken in fast read-out mode. Bias and flat-field 
images were also observed on the same nights, except for the nights 
2008-10-03 and 2008-11-21. 

Each target field in our program was imaged in the g',r', and 
i! filters. For each filter we had three exposures of 180 sec, which 
were slightly dithered by ~ 10" in the direction perpendicular to 
the gap between the SOI CCDs. This dithering pattern allowed us 
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Figure 2. Distribution of main sample clusters as a funtion of photometric 
redshift. 

to fill the gap region as well as to remove CCDs defects and cosmic 
ray hits. The centre of the target cluster was placed at 30" from the 
gap, to the east direction. 

To achieve the image quality needed for the survey in or- 
der to increase the arc detection efficiency, we required a seeing 
FWHM smaller than 0.8". We also required an airmass constraint 
of X 1.5, which is adequate for targets close to the celestial 
equator imaged from SOAR. Finally, we required nights with ^ 7 
days from New Moon, to reduce the noise over the images, mainly 
in the g' band. 



3 DATA REDUCTION 

The individual exposures of each SOGRAS field were bias- 
subtracted and flat-fielded using standard tasks from the MSCRED 
package of the Image Reduction and Analysis Facility (IRAF). For 
the nights in which there were no bias and flat-field exposures, ex- 
posures of the previous night were used. 

Custom codes were employed to remove the fringe pattern 
from i' band images, since the defringing performed by IRAF MS- 
CRED tasks was not satisfactory and we found that it could be 
improved upon. Our codes identify pixels on top of the fringing 
pattern and compute median counts for them, both for the sky sub- 
tracted science and fringing correction images, computing the ratio 
between the two. They then scale the correction frame by this ratio 
before subtracting the pattern from the science image. The fringe 
amplitude (i.e., the typical difference between peak and valley in 
the fringe pattern) in the raw data was approximately 4% of the 
sky and was reduced to levels smaller than 1% of the sky by our 
defringing method for most of our i! band images. However, we 
noticed a remaining fringing residual, especially in the west CCD 
chip, for about 12% of our i' band images. This residual fringing 
amplitude was at worst ~ 2% of the sky level. We also noticed a 
small difference in the counts of the four SOI amplifiers (< 2%) 
and a difference in the noise level between the two CCDs. 

The original Multiple Extension Format (MEF) files were then 
converted into Flexible Image Transport System (FITS) files using 
the task soimosaic from the SOAR/SOI IRAF package. 

Exposures of the same filter of each target were aligned and 
combined into a stacked image, by taking the median at each posi- 
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Figure 1. On-sky distribution of SOGRAS clusters. The area bounded by the dashed line is the CS82 footprint, which essencially is a subarea of the Stripe 
82 footprint (bounded by the solid line). The auxiliary sample is outside of Stripe 82, but is still in the main SDSS footprint. For scale reasons, the 3 clusters 
whose positions are far from the equatorial stripe (SGAS0940+0744, SGAS1023+0413, and SGAS 1054+ 1439) are not shown. 















in 

14 








1 
1 
1 
1 


























12 

I 10 
S 

6 




















1 












1 






2 































ftl i-i- . . L_,_J : ' ! - I 

0.4 0.G 0.8 1.0 1.2 1.4 1.6 



/ I! //.'..' (arcsec) 
Figure 3. Seeing FWHM distribution of the SOGRAS fields. 



tion. This stacked image was used for object detection and photom- 
etry. The stacked images in each filter were aligned with the task 
wregister from IRAF, using the i' band image as reference. We also 
combined the g', r', and i' stacked images of each cluster to have 
a final g' + r' + i' coadded image. This latter was used to visually 
inspect for gravitational arc candidates. 

We measured the seeing in the stacked images for each band 
using the imexamine IRAF task. The distribution is shown in Fig. [3] 
Clearly most of the images do satisfy our seeing constraint, at least 
in the i' band, and are therefore well suited for finding arcs. Only 
two clusters have seeing larger than 1" in that band. They were ob- 
served in windy conditions and their images have relatively poor 
seeing. The median seeing for all images is 0.82", 0.74", and 
0.69" in the g' , r' , and i bands respectively. 



3.1 Astrometric calibration 

For the construction of the world coordinate system (WCS) of each 
stacked image, we used The Guide Star Catalog, Version 2.3.2, 
GSC2.3-STSCI JLasker et al.|2008> and pattern-matched the posi- 
tions of the stars against those in the SOGRAS fields, producing a 
list with the RA, Dec, and cartesian (x and y) coordinates of the se- 
lected stars. For this process we selected a subsample of the bright- 
est (but not heavily saturated) stars yielding typically about 20 stars 
for each SOGRAS field. The residuals in the astrometric solutions 
for RA and Dec are of the order of 0.15", which is sufficient to 
enable a proper positional matching to SDSS for the photometric 
calibration. 

3.2 Photometric calibration 

We used public SDSS data to calibrate the photometry of our 
stacked images. Following the standard SDSS system, we hereafter 
refer to our calibrated data as g, r, and i magnitudes. Most fields 
we observed belong to the SDSS Stripe 82 and therefore have eas- 
ily available and very accurate photometric calibrations ftvezic et] 
|al.|2007) . 

The technique we employed to calibrate our sample was as 
follows. The first step was to select bright but unsaturated stars 
from SDSS DR7 jlvezic et al.] [20071 in the SOGRAS fields. 
For each field, we detected objects with SExtractor Version 2.8.6 
( |Bertin & Amouts|1996} and defined as stars objects whose star- 
galaxy classification parameter CLASS-STAR, as measured by 
SExtractor in the i band image, was larger than 0.85. Those were 
matched to the stars in the SDSS catalog, yielding typically about 
10 stars per SOGRAS field. Then we computed the mean off- 
set between the magnitudes measured with SExtractor's automatic 
aperture (MAG^UTO) and their corresponding SDSS magnitudes 
(MODELJVIAG), and used this offset as the zero-point of the mag- 
nitude scale of each filter for that field. 

Given the lack of actual photometric standards in our calibra- 
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Table 1. Summary of the SOGRAS galaxy cluster sample. 



Cluster ID RA Dec 

z phot Observation date 

(J2000) (J2000) 



a Cluster from the auxiliary sample. 
6 Cluster from the extra sample. 



tion process and the fact that our focus is on extended objects, we 
refrained from adopting colour terms in the calibration, whose am- 
plitude is likely smaller than the uncertainties in the galaxy pho- 
tometry. 



3.3 Galaxy photometry 

We used SExtractor to identify sources and measure their magni- 
tudes. Based on provided detection thresholds, SExtractor deter- 
mines the background around each source and whether a given 
pixel belongs to the source or to the background. SExtractor au- 
tomatic aperture photometry (MAG_AUTO) was adopted in this 
work. It is based on flexible elliptical apertures around every de- 
tected object and is intended to give the a precise estimate of total 
magnitudes, at least for galaxies. 

An automated pipeline in python was created to expedite the 
process of object finding, photometry (including application of 
photometric calibration described earlier), and catalog construc- 
tion. This pipeline uses functions from SLtool^] a library for 
image processing, catalog manipulation, and strong lensing appli- 
cations ( |Brandt et al.| [2012). In the following we provide a brief 
description of this pipeline. 

For each SOGRAS field, our pipeline runs SExtractor sepa- 
rately for the g, r, and i stacked images, using the previously de- 
termined zero-points. Sources with more than 10 pixels and whose 
flux exceeds 2a above the sky background were considered as real 
detections^ 

In section|4]we describe the resulting catalogs in more detail. 

3.4 Photometric Quality Assessment 

In order to assess the quality of the SOGRAS photometry, we again 
used the well-calibrated data from SDSS, given the overlap of the 
SOGRAS fields with the SDSS footprint. We began by performing 
an object matching for each SOGRAS field, as was done in the 
selection of stars for photometric calibration. 

We compare our magnitudes and colours to those from SDSS 
in order to assess our photometric errors and their dependence on 
S/N level. In Fig. [5] we show, as an example, the comparison of 
i band magnitudes and (r — i) colours for both stars (shown on 
the left panels) and galaxies (right panels) in the field of the clus- 
ter SGAS0850+0015. At the bright end, the scatter in the plots is 
dominated by the differences in the way magnitudes were measured 
in SOGRAS and SDSS, and by the residuals in photometric cali- 
bration. On the other hand, at the faint end of the plots, the larger 
scatter is probably caused by the low S/N levels of these objects, 
specially in SDSS, which is shallower than SOGRAS. No signifi- 
cant systematics is seen in the stellar photometry, indicating that the 
photometric calibration is effective. After applying a 2.5-a clipping 
to eliminate outliers, we find a mean offset of {i — isDSs) = 0.02 
for i < 19. The rms SOGRAS-SDSS residual is i = 0.04 in the 
same range. Stellar colours have slightly larger systematic residu- 
als (r — i) — (r — i)sDSS = 0.05 with an rms value of 0.04 when 
bright stars are considered. These values are typical of the other 
fields. 

As for the galaxies, the mean offsets in i magnitudes and 
(r — i) colours depicted in Fig.|4]are comparable to those of stars. 
But the scatter is larger as can be attested from a visual inspection 
of Fig. [4] This reflects the difficulty in measuring their total fluxes, 
specially for faint objects, and also de differences in the SOGRAS 
and SDSS PSFs. The rms residuals are 0.06 both for magnitudes 

6 The SLtools library is available at 

http : // che . cbpf .br/sltools/ 

7 The remaining parameters used for object detection and catalog genera- 
tion with SExtractor can be obtained from the configuration file, which is 
available upon request to the corresponding author. 
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Figure 4. Differences in j band magnitude (upper panels) and (r — i) colour (lower panels) between SOGRAS and SDSS for both stars (left panels) and 
galaxies (right panels) in the field around the cluster SGAS0850+0015. Only objects with r < 22.5 were used in the plots. The solid line corresponds to zero 
residual. The square points were 2.5<r clipped out before computing the mean residual (dashed line). 



and for colours. The scatter plot in colours also reveals a system- 
atic trend of SOGRAS galaxies being redder (bluer) for blue (red) 
colours, which is another way of concluding that the SOGRAS 
colour distributions are narrower than those from SDSS. This is 
just as expected from the higher photometric precision and image 
quality from our SOAR data. 

We now extend our photometric quality assessment to the 
other SOGRAS fields. In Figure [5] the panels on the left show 
the distribution, over all fields, of the mean difference between our 
magnitudes and those from SDSS using bright sources only. The 
panels on the right show the rms difference, computed as follows: 



1 

71 



(rrij - m^sDSs)' 



N 



1/2 



(1) 



where rrij are the individual magnitude measurements from either 
surveys in a given field and N is the number of bright sources 
in that field. By bright sources in this figure we mean those with 
g < 20, r < 20, i < 19. The y/2 factor in the expression for the 
rms accounts for the fact that errors in both SDSS and SOGRAS 
contribute in quadrature to this statistic. Therefore, we are here tak- 
ing the relatively conservative approach of that both photometric 
samples are of equal precision, which may lead to an overestimate 
of the SOGRAS random photometric errors. The solid (dashed) his- 
tograms in each panel are for the stars (galaxies). Each row corre- 



sponds to a given filter. Only fields with at least 5 bright sources 
were included in the histograms. 

The mean photometric offsets for the stars are well centered 
around (m — msDSs) = 0, with very few SOGRAS fields having 
(m — msDSs) | > 0.1. The global median values of these offsets 
are 0.007, 0.007, and 0.005 for g, r and i, respectively. This result 
shows that our calibration was successful and consistent for the 
whole survey. The mean offset for galaxies is clearly larger and 
systematically positive. However, the typical galaxy systematics is 
still constrained to (m — msDSs) Ss 0.2 in most cases. This larger 
systematics reflects the complexity of measuring galaxy fluxes and 
is likely caused by the differences in measuring methods applied to 
SOGRAS (MAG_AUTO) and SDSS (MODEL JVI AG). Notice that 
the mean galaxy offset also tends to be larger in the i band, likely 
as the result of fringing residuals accummulated on galaxy angular 
scales. 

The rms plots bear information on the random rather than 
systematic effects. Since they are also restricted to bright objects, 
we can estimate the photometric uncertainty in our calibration us- 
ing the stellar rms of a typical field: rms(g) ~ 0.03; rms(r) ~ 
0.03; rms(i) ~ 0.03. For the galaxies, these typical values are 
larger, rms(g) ~ 0.11; rms(r) ~ 0.08; rms(i) ~ 0.11, since 
they incorporate the effect of different magnitude definitions on top 
of that from calibration. 

In Figures |6]|7j and [8] we assess the SOGRAS random pho- 
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Figure 5. Panel (a): distribution of mean g band magnitude differences between our photometry and that of SDSS for bright stars (solid) and galaxies 
(dashed). Panel (b): distribution of rms g band magnitude differences between our photometry and that of SDSS for bright stars (solid) and galaxies (dashed). 
See equation^and text for details. Panels (c) and (d): same as in (a) and (b) but now for the r band. Panels (e) and (f): same as in (a) and (b) but now for the i 
band. 



tometric errors as a function of S/N level. The values of Sm are 
computed as 



8m — 



Ns (m) , 



N s (m)V2 ' 



(2) 



where the sum is over all Ns(m) SOGRAS fields with at least 
5 galaxies in that m bin. rmsj is the root mean square residual 



relative to SDSS among galaxies in the j field at that magnitude 
bin. The errorbars in the figures are the dispersion around the Sm 
in each magnitude bin. 

The figures show a systematic increase in the photometric un- 
certainties as fainter galaxies are considered. But notice that a min- 
imum S/N — 3 (which corresponds to Sm = 0.36) is achieved for 
all but the faintest magnitude bin. From these plots we infer that, in 
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Figure 6. Magnitude uncertainties in g estimated according to Eq.[2]shown 
in the text. This is basically an average (over all SOGRAS fields) of the 
rms residuals relative to SDSS divided by v2 in each bin. The error bars 
are the dispersion around this average. The 8 m = 0.36 line indicates the 
S/N = 3 level. 
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Figure 7. Same as Figure|6] but for the r band. 

average, higher S/N are achieved in the g, than in the r, and then 
in the i bands, with detection limits for S/N > 3 being roughly 
g~23, r~22.5, and i ~ 22. 



4 GALAXY AND ARC CATALOGS 
4.1 Galaxy catalogs 

We merged the catalogs resulting from the photometry (see sec- 
tion [3T3j of the 3 bands in a final object catalog, choosing the i 
band catalog as reference for object position. This means that we 
searched in the g and r band catalogs for the objects that are in the 
i catalog. Non-matched objects received a flag —99.99 in the cor- 
responding filter. The final catalog was saved as a FITS table and 
contains the information on position, magnitudes, morphological 
parameters, and star-galaxy classification of the detected objects. 
In Table [2] we list the parameters in the final object catalog and 
their corresponding definitions. This catalog is publicly available 
upon request to the corresponding author. 



Figure 8. Same as Figure[6] but for the i band. 



Table 2. Column labels in the SOGRAS object catalog. 



Parameter name 


Definition 


CLUSTERJD 


Cluster ID 


OBJECT JD 


Object number 


RA a 


Right ascension (J2000) 


DEC 


Declination (J2000) 


X_IMAGE a 


Object position along x 


YJMAGE" 


Object position along y 


MAG_AUTO_G 


Automatic aperture magnitude in the g band 


MAGERR^JJTCLG 


RMS error for AUTO magnitude in the g band 


MAG_AUTO_R 


Automatic aperture magnitude in the r band 


MAGERR^JJTO-R 


RMS error for AUTO magnitude in the r band 


MAG_AUTO_I 


Automatic aperture magnitude in the i band 


MAGERR^UJTOJ 


RMS error for AUTO magnitude in the 2 band 


THETA.SKY 6 


Position angle (east of north) 


AJMAGE 6 


semi-major axis for second moments 


BJMAGE 6 


semi-minor axis for second moments 


ELLIPTICITY 6 


Ellipticity (1 - B JMAGE/AJMAGE) 


CLASS _STAR a 


Star-galaxy classification 


FLAGS" 


Extraction flags 


"Measured in i band 



Measured in r band 



In Figure|9]we show colour-magnitude diagrams for two clus- 
ters, one at the low-z bin (SGAS0850+0015, z = 0.20, upper pan- 
els) and the other with higher z (SGAS0202-0055, z = 0.50, lower 
panels). The red sequence is more visible in the lower z cluster, 
as expected. The colours of the red sequence also change slightly 
with z. The variation seems to be systematically larger for (g — r). 
For the relatively few high-z clusters where the red sequence is 
clearly visible, it tends to have bluer colours (g — r ~ 0.6) than 
the low-2 ones, for which (g — r) > 1. This is likely caused by the 
4000 A B aimer break affecting the g filter at z < 0.4 but not the 
other passabands. 

We chose to measure the morphology — semi-major (A) and 
semi-minor (B) axes and position angles derived from weighted 
second moments in SExtractor — in the r stacked images to op- 
timize the combination of S/N and seeing. The seeing is not sig- 
nificantly degraded from the i to the r band (see Fig.[3j, while the 
S/N increases from most objects (c.f. section [3!4|. Therefore, the 
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Figure 9. Colour-magnitude diagrams of clusters SGAS0850+0015 (upper panels) and SGAS0202-0055 (lower panels), (g — r) colours are shown on the left 
panels whereas (r — i) colours are on the right panels. The triangles indicate galaxies that are located in a circular region of 100" (upper panels) and 60" 
(lower panels) around the cluster centre. 
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Figure 10. Left panel: ellipticity distribution of all SOGRAS galaxies (all galaxies) and of SOGRAS galaxies which are located close to the clusters centre 
(cluster galaxies). Right panel: ellipticity distribution of SOGRAS galaxies which are located close to the clusters centre for low and high z clusters. Galaxies 
were selected as objects with star-galaxy parameter CLASS.STAR < 0.85, as measured in r band images. 
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r band provides a balance between the higher S/N of the g band 
and the better seeing of the i band. Figure [10] shows the distribu- 
tion of the ellipticity e = 1 - B/A of the SOGRAS galaxies. The 
left panel shows the distribution of e for all objects classified as 
galaxies and for all of them which are located close to the centre of 
the SOGRAS clusters. Both distributions are peaked at e ~ 0.15, 
and are strongly skewed towards higher values. This behaviour is 
qualitatively confirmed by analyses of early-type galaxies in nearby 
clusters (e.g., Fasano et al. 20101. In the right panel we present the 
distribution of e for SOGRAS galaxies that are located close to the 
clusters centre for low and high z clusters. The similarity between 
the distributions for low and high z clusters indicates that our shape 
determinations for high-z objects are not strongly degraded by at- 
mospheric seeing and that shapes were properly measured for the 
galaxies in our catalog. 

4.2 New arc system candidates 

We inspected all SOGRAS images in order to look for strong 
gravitational lens systems. We found 6 clusters (SGAS032 1+0026, 
SGAS0328+0044, SGAS00 14-0057, SGAS0041-0043, 
SGAS0940+0744, and SGAS1023+0413) that show clear ev- 
idence of arcs and 2 clusters (SGAS0219+0022, SGAS0202-0055) 
that show probable arcs. 

We identified 16 arc candidates close to the brightest mem- 
bers of 8 cluster cores and most of them show bluer colours than 
the central cluster galaxies. Four of them are giant arcs and have 
length- to- width ratio (L/W) larger than 7. The remaining candi- 
dates are arclets, i.e. have smaller L/W. 

The position, magnitude, length, and width of the arc candi- 
dates in the 6 most probable strong lenses are displayed in Table 
[3] Contrary to the measurements of magnitude of other objects in 
the catalog, which were obtained with SExtractor, the magnitudes 
of the arc candidates were measured using the task polyphot from 
IRAF. This task computes the magnitudes inside polygonal aper- 
tures, providing more precise measurements of arc magnitudes, 
since it takes into account the arc shape. The length and width of 
the arc candidates were visually estimated. 

In Figure [TT] we show the candidate strong gravitational lens 
systems identified in SOGRAS data. Notice the systematically 
bluer colours of the arcs in comparison to that of a typical cen- 
tral cluster galaxy in our range of redshifts, (g — r) ~ 0.8 — 1.0 
(see Tableland online colour version of Figure[TTJ. 

From the 8 lens system candidates, 2 of them are in the low-z 
bin, 4 of them are in the high- z bin and the remaining two are in 
the extra sampl^] These later cannot be used in statistical analy- 
sis, since the selection criteria for the extra sample are not the same 
as for the low and high z samples. Thus, the efficiency of arc for- 
mation inferred from SOGRAS is 4/25 = 16% ± 8% for high-z 
clusters and 2/22 = 9% ± 6% for low-z ones, including only Pois- 
son errors. 

Concentrating on the 6 lensing systems from the low and high 
z samples, we infer that about 10% of the clusters have arcs around 
them. This overall efficiency is in agreement with larger arc sur- 
veys, such as |Gladders et al.| ^2003} and |Hennawi et al.| ( |2008| >. 
Despite the low number statistics, the results are also in qualita- 
tive agreement with the models in that they predict the high-z bin 

8 This result indicates that the clusters with x-ray match are more efficient 
lenses than the overall clusters selected for the main SOGRAS sample, as 
expected (e.g.JHoresh et al.|201 11. 



to have a larger efficiency in arc formation (e.g., jCaminha et7aT| 
|20"l2l >. 

Follow-up observations of the first 3 lens systems candi- 
dates identified in SOGRAS images (in clusters SGAS0321+0026, 
SGAS0328+0044, and SGAS02 19+0022) were conducted on the 
8 m Gemini Telescope with the Gemini Multi-Object Spectrograph 
( Cypri ano et al.|20T0| >. The main aims of this follow-up programme 
were to confirm spectroscopically the gravitational lensing nature 
of these candidates, provide mass estimates for the clusters from the 
velocity dispersion of their member galaxies, and perform strong 
lensing reconstruction of the projected mass distribution of the 
lenses. The results will be shown and discussed in a forthcoming 
paper. 

Besides the multi-object spectroscopy, we obtained deep 
imaging to search for new arcs and to determine properties of 
the lensed galaxies (sources) such as their stellar populations and 
star formation rate. Visual inspection on these deeper images con- 
firms the arc candidates found in clusters SGAS032 1+0026 and 
SGAS0328+0044, but revealed that the candidates of the cluster 
SGAS02 19+0022 are unlikely to be arcs. The two of them which 
lie closer to the cluster are clearly seen as point sources in these 
images, while the third revealed itself as a superposition of two 
relatively edge-on galaxies. On the other hand, from a visual in- 
spection of these images, several arc candidates are found in all of 
them, including in SGAS0219+0022. 

We thus conclude from our preliminary analysis of the deeper 
Gemini data that the number of lensing clusters has not been 
changed relative to our SOAR based search, although some individ- 
ual arc candidates have been added and some removed. Therefore, 
our preliminary analysis of arc statistics above is kept unaltered. 



5 SUMMARY AND FUTURE PERSPECTIVES 

We presented the first results from SOGRAS, an imaging survey to- 
wards 5 1 galaxy clusters using the SOAR telescope. We carefully 
assessed the quality of our data products. We estimate our galaxy 
detection limits as g ~ 23, r ~ 22.5, and i ~ 22 at S/N ~ 3. 
Photometric calibration was performed using the SDSS stars in 
common with our fields with systematic uncertainties amounting 
to 0.007, 0.007 and 0.005 mag in g, r, and i respectively, and 
with essentially no systematics. Galaxy photometry suffered from 
a systematic offset typically of 0.04 — 0.12 mag in comparison to 
SDSS, likely caused by the different measurement methods used. 
Our source catalog has over 21,000 entries, about 90% of which are 
galaxies. We confirm that the data are of enough quality to allow a 
clear red sequence to be seen in most clusters in the low-z bin. Fur- 
thermore, seeing effects have not strongly affected shape measure- 
ments from our images, as attested by the fact that the distribution 
of axis ratios from our data closely resembles that for nearby cluster 
galaxies. 

Although the number of targets in SOGRAS is smaller than 
in other arc surveys, the strength of our survey resides in the fo- 
cus on two narrow redshift intervals whose differential lensing ef- 
ficiency may yield direct information about the evolution of arc 
incidence. This works in a complementary way to previous studies 
of arcs around galaxy clusters, some of which have a larger overall 
statistics in the complete survey, but not in these redshift bins (e.g., 
|Hennawi~et al. 2008 1. The observations where carried out in similar 
conditions, with very good seeing, and with the same instrument, 
assuring an homogeneity of the data. Furthermore, our sample is 
unbiased, in the sense that there was no selection based on an a 




Figure 11. (Color online) Strong lensing candidates identified in SOGRAS data. In all cases, the left panel is the g + r + i coadded image and the right 
panel is the colour-composite image. 
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Figure 11. — Continued 
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Table 3. Properties of the arc candidates identified in SOGRAS data. The photometric redshift z p ^ ot of the clusters was taken from Table 1. 



Cluster ID 


Zphot 


Arc ID 


RA 
(J2000) 


Dec 
(J2000) 


9 


V 


i 


L 

(arcsec) 


W 
(arcsec) 


L/W 


SGAS0321+0026 


0.309 


A 


03:21:10.55 


00:26:20.64 


22.84 


22.41 


22.61 


5.08 


1.08 


4.71 


B 


03:21:12.76 


00:26:23.73 


24.13 


23.63 


22.90 


1.23 


0.62 


2.00 


SGAS0328+0044 


0.322 


A 


03:28:15.79 


00:44:49.59 


22.73 


20.72 


20.25 


5.93 


3.03 


1.95 


B 


03:28:15.34 


00:44:57.12 


24.29 


21.84 


21.39 


2.70 


1.39 


1.94 


SGAS0014-0057 


0.535 


A 


00:14:54.93 


-00:57:02.44 


21.84 


21.21 


20.65 


12.24 


1.31 


9.34 


SGAS004 1-0043 


0.564 


A 


00:41:09.22 


-00:43:47.38 


22.65 


21.57 


20.89 


9.16 


1.40 


6.54 


SGAS0940+0744 


0.390 


A 


09:40:53.33 


07:44:17.48 


23.52 


21.53 


21.95 


11.70 


1.23 


9.51 






A 


10:23:38.59 


04:11:20.77 


24.43 


23.17 


22.41 


2.77 


0.8 


3.46 


SGAS1023+0413 


0.465 


B 


10:23:41.45 


04:10:41.62 


24.36 


23.19 


22.57 


2.93 


1.10 


2.66 






C 


10:23:41.43 


04:11:00.11 


23.91 


23.44 


22.71 


2.39 


0.93 


2.57 



priori likelihood for an individual cluster to have arcs. These fac- 
tors make the SOGRAS sample well suited for arc statistics studies, 
despite the relatively low number of objects. 

Preliminary results from a visual inspection suggest an overall 
efficiency for arcs of about 10%, consistent with previous studies 
pennawi et al.|20 08 ; Gladd ers et al.|2003| >. Furthermore we find 
a possible indication of a higher arc incidence in the high-z bin 
with respect to the low-z, qualitatively in agreement with theoret- 
ical modeling (Caminha et al. 2012). A detailed study of arcs in 
SOGRAS, including comparison with model predictions, mass es- 
timates from arcs, quantitative studies on arc morphology and arc 
detection will be presented in a separate paper. 

Besides the strong lensing studies, this good image quality 
data can be used to perform a high signal-to-noise weak-lensing 
analysis by stacking the weak lensing signal of all clusters in a 
given redshift bin to obtain an overall mass estimate for the clus- 
ters. We will also use the arcs and other strong lensing features to 
constrain the individual masses of the clusters (e.g., Cypria no et aL] 
|2005[ >. The data will also be used for detailed galaxy morphological 
studies using model-fitting methods, including the modeling of the 
PSF. 

As far as we know, this is the first arc survey that specifically 
targeted optically selected clusters from the deep Stripe 82 coadd. 
This enabled us to select our low-z and high-z samples. Further- 
more, this opens the possibility to exploit the combination with 
other surveys, by cross matching with the wealth of data in that 
region of the sky. In particular, most fields are in the CS82 region, 
which will provide complementary information around the clusters 
on larger scales. 

This dataset will also be used to validate and benchmark arc 
identification and characterization tools being developed by our 
group, including methods to enhance their detectability and arc 
finding algorithms. In some aspects, this small survey can be seen 
as a pathfinder for stage III photometric surveys such as DES, given 
that the SOAR images have similar depth as expected for those sur- 
veys. 

SOGRAS can also be seen as a feasibility study for a unique 
arc survey profiting from the adaptive optics capabilities of SOAR 
in the optical and NIR, which can improve the PSF by a factor of 2 
to 5. 
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